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Analytical  Characterization  of  Bistatic  Scattering 
From  Gaussian  Distributed  Surfaces 


1.  INTRODUCTION 


Accurate  statistical  modeling  of  the  clutter  encountered  by  a  radar  system  at  a  particular 
site  Is  critical  to  the  proper  design  of  the  system.  Clutter  Is  any  unwanted  radar  echo,  such  as 
reflections  from  the  terrain.  When  the  clutter  Is  as  large  as  the  expected  target  return,  the  radar 
assumes  the  return  Is  due  to  a  target  being  present,  and  a  false  alarm  occurs.  Clutter  Is, 
therefore,  a  significant  factor  In  determining  the  optimum  value  of  the  minimum  detectable 
signal,  (smallest  target)  for  the  required  false  alarm  rate.^  The  clutter  returns  are  not  equal  but 
are  distributed  in  magnitude.  If  there  are  more  high  energy  returns  due  to  clutter  than  a  model 
predicts,  the  false  alarm  rate  will  be  unacceptably  high.  Conversely,  fewer  high  energy  returns 
than  expected  means  that  the  probability  of  detection  could  have  been  improved  by  lowering 
the  detection  threshold. 

Also,  theoretical  modeling  of  scattering  from  rough  surfaces  that  cause  clutter  Is  very 
Important  because  accurate  measurements  of  clutter  are  very  difficult  and  expensive  to 
perform,  and  are  sometimes  unreliable.  The  mean  and  variance  of  clutter  can  be  determined 
experimentally,  but  this  requires  many  measurements  of  similar.  Independent  clutter  cells  to 
be  made.  These  measured  values  would  then  only  be  vabd  for  the  specific  geometry  and  type  of 

Received  for  Publication  19  November  1991 

'  Skolnlk.  M.l.  (1980)  Introduction  to  Radar  Systems,  McGraw-Hill  Book  Company,  New 
York. 
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clutter  being  measured.  Any  changes  In  polarization,  frequency,  radar  deployment,  geometry, 
seasonal  variation,  or  even  moisture  content  of  the  ground  will  result  In  a  completely  different 
set  of  statistics.  Even  measuremenls  made  under  Identical  conditions  can  vary  by  up  to  10  dB.^ 
For  these  reasons,  theoretical  models  of  clutter  statistics  are  Important.  The  power  scattered 
from  a  rough  surface  Is  usually  assumed  to  be  Rayleigh  distributed.  The  Rayleigh  distribution 
function  Is; 


f(x) 


2  2 

—  e^*“  U(x). 
a 


(1) 


Figure  1  Is  a  plot  of  this  distribution  function  for  several  values  of  a.  The  assumption  of 
Rayleigh  distributed  clutter  Is  valid  for  an  Infinite  scattering  surface  of  many  similar  sized 
scatterers.  If  the  scatterers  differ  significantly  In  size  or  If  the  number  of  scatterers  Is 
decreased  by  limiting  the  cell  size,  the  probability  density  function  of  the  scattered  power  will 
no  longer  have  a  Rayleigh  distribution. 

An  application  of  this  Is  In  short  pulse  radars  where  the  clutter  patch  is  limited  In  the 
range  dimension.  For  such  radars,  the  requirement  of  a  laige  surface  with  many  scattering 
centers  Is  violated.  In  this  repiort,  the  statistics  of  the  power  scattered  from  a  rough  surface 
with  a  Qnlte  range  dimension  will  be  determined  using  a  physical  optics  model. 

The  effect  of  polarization  on  rough  surface  scattering  will  also  be  examined.  Because  the 
reflection  coefficient  of  a  surface  Is  polarization  dependent,  the  statistics  of  the  scattering  wlU 
also  be  a  function  of  the  Incident  polarization.  In  a  study  by  Papa  and  Woodworth,^  the  mean 
and  variance  of  the  scattered  power  were  Investigated  for  horizontal  Incident  polarization. 
Arbitrary  linear  polarization  states  will  be  studied  to  determine  the  effect  of  polarization  on 
the  scattering  statistics. 


^  Long,  M.W.  (1975)  Radar  Reflectivity  of  Land  and  Sea,  D.C.  Heath  and  Company, 
Lexington,  MA. 

^  Papa.  R.J.,  and  Woodworth.  MB.  (1991)  The  Mean  and  Variance  of  Diffuse  Scattered 
Power  as  a  Function  of  Clutter  Resolution  Cell  Size,  RADC-TR-91-09. 
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Figure  1.  Rayleigh  Distribution  Fi 


2.  THEORETICAL  FORMULATION 


The  normalized  cross  section  of  a  rough  surface.  o°  is 

['  ( VE,)  >  -  <  (V^i) ]  12) 


where  E,  Is  the  Incident  field.  Eg  Is  the  scattered  field,  Rq  Is  the  distance  Ifom  the  surface  to  the 
field  point,  A  Is  the  area  of  the  Illuminated  cell,  *  indicates  the  complex  conjugate,  and  the 
angle  brackets  Indicate  an  ensemble  average  over  the  surface  variables.  Figure  2  depicts  the 
scattering  geometry.  For  a  one  dimensionally  rough  surface  only  in-plane  scattering  occurs 
|<t»  =  0°),  and  there  Is  no  shadowing.'^ 


Fljiure  2.  Scattering  Geometry 


*  Beckmann.  P..  and  Splzzlchlno,  A.  (1963)  The  Scattering  of  Electromagnetic  Waves  from 
Rough  Surfaces.  The  MacMillan  Co.,  New  York. 
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The  surface  geometiy  is  shown  In  Figure  3.  The  points  xj  and  X2  are  two  random  points 
separated  by  the  correlation  distance,  t.  with  heights  and  ^2-  The  correlation  distance 
defines  the  density  of  the  Irregularities.  The  local  slopes  of  the  rough  surface  at  these  points 
are  g|  and  H2-  dlsliibutlon  of  the  heights  is  described  by  a  Gaussian  distribution  function 
liecause  this  is  the  most  typlc  id  distribution  of  a  rough  surface.  This  distribution  function  Is 
given  by: 


w  (^) 


-J—  f-xp 
a  V2x 


2  o'  J 


(3) 


where  o  Is  the  standard  delation  of  the  surface  heights  and  o^  Is  the  variance  of  The 
variance  can  be  chosen  to  represent  the  desired  degree  of  roughness  and  the  correlation 
distance  chosen  for  the  desired  density  of  the  Irregularities. 

For  a  one  dimensionally  rough  surface  used  In  this  study,  the  expression  for 
normalized  mean  cross  section,  o”.  Is: 


o 

o  = 


(4) 


The  scattered  field  Is  a  cylindrical  wave  calculated  by  evaluating  the  Helmlioltz  Integral 
equation: 


E 

S 


V'i?.  dL 
9n, 


(5) 


where:  V  =  V  2/()Tk^  •  R ')  exp  (1  •  R') 

ks  =  scattei  ed  wave  vector 
R  =  Rq  —  Tq 

Ro  =  vector  from  origin  to  field  point 
r  =  vector  from  origin  to  source  on  the  surface 
E  =  (1+R)  E|  Is  the  total  field  on  the  surface 
R  =  Fresnel  reflection  coefficient 

--  =  normal  derivative  on  surface,  S 
on 

L  =  ramge  dimension  of  the  scattering  cell 
and  the  time  dependence,  exp  (-Jcot).  Is  supressed. 
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Figure  3.  Surface  Geometry 


The  most  general  expression  for  o"  Involves  a  six- fold  Integral  over  the  surface  variables, 
xi,  X2.  ^1.  ^2-  However,  Papa  and  Woodworlfi  show  that  because  the  surface  heights 

are  regarded  as  a  stationary  random  process,  the  Integral  Is  a  function  solely  of  the  separation 
between  the  two  surface  points,  t  =  xi  -  X2,  which  simplifies  the  six-fold  Integral  to  a  five-fold 
Integral.  The  flve-fold  integral  can  be  further  simplified  to  a  four-fold  Integral  because  the 
expression  for  o*‘  is  a  function  of  the  height  differences,  4  =  ^i  ~  42-  Then,  by  assuming  that  the 
surface  slopes  are  small.  V2  o/  T  <  1.  the  Integrals  over  pi  and  P2  can  be  performed  analytically 
and  the  expression  for  n°  cam  tlien  be  written  as  a  double  integral  given  by: 


o 

0  = 


^  I 

2LXi 


pL/2 

dx^ 

1 

./-L/2  . 

/-L/2 

dx^^Xa-XjXi.cos  (v^  x) 


(6) 


where:  x  =  xi  -  X2 

X2  =  exp  (-l^  (1  _  Ci2)j  Is  the  chairacterlstlc  function  for  bivariate  height  distribution 
Xi  =  exp  (-£^/2)  Is  the  characteristic  function  for  univariate  height  distribution 
Ci2  =  exp  Hxi  -  X2)2/'I^1 

V  ^  =  Rayleigh  parameter  squared 
o  =  rms  surface  helglit 
Vx  =  (2  n/  X)  Isln  lOi)  -  slnOa)) 

V2=  (2  jt/  X)  (cos  (0i)  +  cos(0a)) 

F  = 

cos  0^  +  cos  0^ 


Note  that  It  Is  the  reflection  coefficient  In  the  factor  F  that  causes  o°  to  be  a  function  of  the 
polarization  of  the  incident  wave  and  the  dielectric  properties  of  the  surface. 

The  variance  of  the  scattered  power,  Is  given  by  (see  Papa  and  Woodworth): 


L/2 


0  = 


nF 

V2LX 


IJJJ 


d  Xj  d  x^  d  X3  d  x^  cos  [v^  (Xj  -  x^  +  X3  -  x^)]  .  [x^  -  x,2  X34] 


(7) 


L/2 


where:  X12  =  ^  d  -  ^12)! 

X34  =  exp  (-12  (1  -  €34)1 
C,j  =  exp  (-X|  -  Xj)/T2| 
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For  a  Rayleigh  distribution,  the  variance  Is  equal  to  the  mean  squared.  This  definition, 
=  (o°  Is  the  definition  used  to  cailculate  the  Rayleigh  variance  In  this  study.'^ 


3.  VERTICAL  POLARIZATION 


The  analysis  done  by  Papa,  et  al  Is  limited  to  the  case  of  horizontal  Incident  polarization.  To 
extend  the  analysis  to  Include  vertical  and  arbitrary  linear  Incident  polarization,  new 
reflection  coefliclents  must  be  calculated. 

For  a  perfect  conductor  the  reflection  coefficient  Is  Rh  =  -1  for  horizontal  polarization, 
(electric  vector  perpendicular  to  the  plane  of  Incidence)  and  Ry  =  +1  for  vertical  polarization, 
(electric  vector  In  the  plane  of  Incidence).  Because  the  reflection  coefficient  is  squared  In 
calculating  the  scattered  power,  the  results  for  horizontal  polarization  and  vertical 
polarization  are  Identical  for  a  perfect  conductor  and  will  not  be  repeated  here. 

Many  surfaces,  however,  have  a  complex  dielectric  constant,  and  will  therefore  have 
different  reflection  coefficients  for  vertical  and  horizontal  polarization.  A  complex  dielectric 
constant  Is  expressed  as  e  =  e'+  je".  where  e'  Is  the  capaclvlty  of  the  medium  and  e"  Is  the 
dielectric  loss  factor.  For  a  good  dielectric  the  term  e’  will  remain  fairly  constant  over  all 
radio  frequencies  and  e"  will  be  small.  The  relative  dielectric  constant  Is  defined  as  Er  =  e/Co 
where  £„  Is  the  capaclvlty  of  vacuum.  Table  1  shows  the  relative  complex  dielectric  constants 
for  the  surfaces  which  are  used  In  this  study  as  well  as  the  dielectric  constant  for  glass,  which 
Is  considered  to  be  a  good  dielectric. 


Table  1.  Suriace  Dielectric  Constants 


Surface 

e' 

e" 

Moist  Ground 

30 

2 

Sea  W.iter 

80 

9 

Dry  Sandy  Loam 

2 

1.6 

Glass 

5.25 

0.0115 
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The  Fresnel  reflection  coefficient  for  horizontal  polarization  is  given  by:® 


= 


cos  0 


-Ve  - 


sin  9 


cos  0 


+  Ve  - 


sin  B 


(8) 


where  0  is  the  local  angle  of  incidence.  The  magnitude  and  phase  of  Ru  are  shown  in 
Figure  4.  Results  for  this  case  given  by  Papa,  et  al  will  be  used  for  comparison  with 
results  obtained  for  vertical  polarization. 

For  vertical  polarization  the  reflection  coefficient  Is: 


R 

V 


e  cos  0  -  A  e  -  sin  0 

r _  r _ 


e  cos  9  +  A  e  -  sin  0 

r  r 


(9) 


Figure  5  shows  the  magnitude  and  phase  of  Rv  for  the  dielectric  constants  used. 

In  comparing  the  horizontal  and  vertical  reflection  coefficients.  Important  differences  can 
be  seen.  The  magnitude  of  the  verbcal  reflection  coefficient  has  a  minimum  value  for  a 
particular  angle  of  incidence.  This  angle  is  called  the  Brewster  angle  and  Is  defined  as 
9b  =  tan~^  V^.  Also,  the  magnitude  of  the  vertical  coefficient  is  always  less  than  that  of  the 
horizontal,  except  at  normal  Incidence  where  they  are  equal.  The  differences  In  the  reflection 
coefficients  will  have  a  direct  effect  on  the  scattering  from  a  surface. 


®  Ruck,  G.T.,  Barrlck,  D.E..  Stuart,  W.D..  and  Krlchbaum,  C.K.  (1970)  Radar  Cross  Section 
Handbook,  Vol.  2,  Plenum  Press,  New  York. 
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Phases  of  Horizontal  Reflection  Coefficients 
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Figure  5.  a.  Magnitudes  of  Vertical  Reflection  CoefBclents 
b.  Phases  of  Vertical  Reflection  CoefBclents 


Many  of  the  other  results,  though,  are  valid  for  both  horizontal  and  vertical  polarization. 
Some  of  the  other  aspects  that  are  considered  are  the  statistical  dependence  of  the  scattering 
on  the  angle  of  Incidence  0i,  the  wavelength  X,  the  length  of  the  clutter  cell  L,  and  the  standard 
deviation  In  the  surface  heights  o.  The  surface  roughness  Is  defined  as  o/T  and  V2a/T  Is  the 
rms  surface  slope.  T  Is  the  correlation  length  and  Is  set  to  one  meter  for  this  study. 

There  are  two  factors  that  cause  the  variance  of  the  scattered  power  to  deviate  from 
Rayleigh;  the  number  of  scatterers  and  the  roughness  of  the  surface.  The  number  of  scatterers 
Is  a  function  of  cell  size.  Incident  angle,  and  the  slopes  of  the  rough  surface.  The  roughness  of 
the  surface  depends  on  the  frequency  of  the  Incident  wave  as  well  as  the  surface  parameters. 
How  these  parameters  affect  the  statistics  of  the  scattered  field  will  be  examined. 

One  condition  necessary  for  Rayleigh  scattering  Is  that  the  illuminated  area,  or  clutter  cell, 
be  very  large  with  respect  to  the  correlation  distance,  so  that  It  will  contain  a  large  number  of 
Independent,  nearly  equal,  homogeneously  distributed  scatterers.  Decreasing  the  size  of  the 
clutter  cell  will  decrease  the  number  of  scatterers  and  should  then  cause  the  distribution  of 
to  deviate  from  Rayleigh.  When  there  are  only  a  few  Independent  scatterers  In  a  particular  cell, 
their  relative  positions  or  the  presence  or  absence  of  one  of  the  scatterers  causes  a  significant 
change  In  the  reflected  signal  from  one  sample  of  the  surface  to  another,  resulting  in  a  large 
variance. 

In  Figures  6  through  10,  several  cases  are  shown  and  the  effect  of  changing  the  cell  size  Is 
examined.  The  scattered  power  Is  normalized  to  both  the  Incident  field  and  the  Illuminated 
area,  and  Is  plotted  as  a  function  of  the  scattering  angle.  Each  plot  shows  three  curves.  The 
curve  consisting  of  long  dashes  represents  the  normalized  mean  of  the  scattered  power,  which 
Is  equivalent  to  the  mean  cross  section.  The  normalized  variance  of  the  scattered  power  of  a 
Rayleigh  distributed  scattering  process  Is  the  curve  of  small  dashes.  The  third  curve  with  the 
crosses  Is  the  calculated  normalized  variance  of  the  scattered  power  for  the  parameters  shown. 

Figure  6  Is  for  0i  =  30°.  X  =  0.25  m.  o  =  0.2  m,  Cr  =30  +  J2  with  L  Increasing  from  L  =  3T 
(Figure  6a)  to  L  =  6T  (Figure  6b)  and  finally  L  =  12T  (Figure  6c).  These  parameters  correspond  to 
a  relatively  smooth  patch  of  moist  ground.  Figure  7  Is  similar  to  Figure  6  but  represents  a 
rougher  surface.  Figures  8  and  9  are  for  a  rough  surface  with  a  75°  angle  of  Incidence,  for 
vertical  polarization  and  horizontal  polarization,  respectively.  The  effect  of  the  different 
reflection  coefficients  Is  noticeable  when  comparing  these  two  figures.  i)artlcularly  at  the 
location  of  the  Brewster  angle,  0b  =85°.  In  Figure  10,  0i  =75°.  X  =  0. 1  m  ajid  o  =  0.2  m 
corresponding  to  a  relatively  smooth  surface.  Figure  1 1  represents  a  fall  ly  rough  water  surface, 
tr  =  80  +  J9. 
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SCATTER  ANGLE  (deg) 

Figure  9.  Mean  and  Variance  of  Scattered  Power  for  H-pol  for. 

a.  Oj  =  75".  X  =  0.25  m.  o  =  0,5.  L  =  3T,  £  =  30  +  J2 

b.  01  =  75".  X  =  0.25  m.  o  =  0.5.  L  =  6T,  e  =  30  +  j2 

c.  01  =  75".  X  =  0.25  m.  o  =  0.5.  L  =  12T,  E  =  30  +  j2 


Mm)  3  00  I  j-  I  L{m)  6.00 
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Figure  10.  Mean  and  Variance  of  Scattered  Power  for: 

a.  01  =  75°,  X  =  0.1  m,  o  =  0.2,  L  =  3T.  e  =  30  +  j2 

b.  01  =  75°  X  =  0.1  m.  o  =  0.2,  L  =  6T,  e  =  30  +  J2 

c.  01  =  75°,  X  =  0.1  m,  c  =  0.2,  L  =  12T.  e  =  30  +  j2 


L(mi  3  00 
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Figure  1 1 .  Mean  and  Variance  of  Scattered  Power  for; 

a.  01  =  75°.  ^=0.1  m,  a  =  0.5.  L  =  3T.  e  =  80  +  J9 

b.  01  =  75°.  A.  =  0. 1  m.  o  =  0.5.  L  =  6T.  e  =  80  +  J9 

c.  01  =  75°.  A  =  0. 1  m.  o  =  0.5.  L  =  12T.  e  =  80  +  j9 


In  all  cases,  it  is  seen  that  for  the  smaller  cell  size.  L  =  3T.  the  scattering  is  more  non- 
Rayleigh  than  lor  the  larger  clutter  cells.  L  =  6T  and  L  =  12T.  in  the  scattering  regions  where 
diflferences  in  the  variance  can  be  seen.  The  variance  of  the  scattered  power  is  consistently 
higher  for  the  smaller  clutter  cells  than  it  Is  for  a  surface  that  gives  rise  to  Rayleigh  scattering. 

Figures  12  througli  15  show  a  comparison  of  the  variance  for  the  different  cell  sizes 
normalized  to  the  Rayleigh  value  so  tliat  a  value  of  one  means  the  calculated  value  of  the 
variance  is  equal  to  the  Rayleigh  value.  The  higher  variance  causes  difficulties  In  radar  design 
because  for  a  constant  mean  value  of  clutter  power,  a  larger  variance  corresponds  to  more 
frequent  large  clutter  returns.  This  Is  a  significant  result  that  must  be  considered  In  the  design 
of  a  short  pulse  radau"  system. 

A  noticeable  trend  .een  In  all  of  tliese  plots  Is  that  the  deviation  from  Rayleigh  scattering 
becomes  more  severe  as  0g  approaches  the  backscatter  region,  particularly  when  the  angle  of 
incidence  is  75°.  Since  scattering  facets  must  have  steep  slopes  to  cause  backscatter,  especially 
for  01  =75°.  and  the  distribution  of  the  slopes  Is  the  derivative  of  the  height  distribution,  and 
therefore  also  Gaussian.®  there  will  be  relatively  few  facets  oriented  properly  to  cause 
backscatter.  Once  again,  the  effect  of  having  a  small  number  of  scatterers  Is  to  Increase  the 
variance  of  the  scattered  power,  causing  the  scattering  to  become  more  non-Rayleigh . 

The  second  factor  that  affects  the  distribution  of  the  scattered  power  Is  the  surface 
roughness.  A  wave  that  is  incident  on  a  rough  surface  will  scatter  In  many  directions  and  It  Is 
the  relative  heights  and  slopes  of  the  scattering  surfaces  that  will  determine  the  angular 
distribution. 

For  a  very  smooth  surface,  the  incident  wave  will  be  almost  entirely  reflected  In  the 
specular  direction  defined  by  Snell's  Law.  As  the  surface  roughness  Increases,  the  scattering 
will  become  more  diffusely  distributed  over  the  entire  angular  space.^  Figure  16a  shows  the 
scattering  for  a  very  smooth  surface,  a/T  =  0.05.  It  Is  quite  clear  that  the  majority  of  the 
scattered  power  is  in  a  small  region  around  the  specular  direction.  In  Figure  16b,  o  =  0.2  m 
corresponding  to  a  relatively  smooth  surface,  and  the  peak  of  the  scattering  In  the  specular 
direction  Is  no  longer  well  defined.  Figure  16c  Is  for  a  rougher  surface  with  o  =  0.5  m.  In  this 
case  the  scattering  is  quite  uniform  over  the  entire  region  with  no  dlscemable  peak  In  the 
specular  region. 


®  Papoulls.  A.  (1965)  Probability.  Random  Variables,  and  Stochastic  Processes,  McGraw-Hill 
Book  Company,  New  York. 
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SCATTER  ANGLE  (deg) 

Figure  12.  Comparison  of  Variance  for  Different  Cell  Sizes  for; 
6(  =  30'=.  X  =  0.25  m.  a  =  0.2.  e=  30  +j2 
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Figure  16.  Mean  and  Variance  of  Scattered  Power  for: 

a.  01  =  30°.  X.  =  0.1  m,  o  =  0.05,  L  =  3T.  e=  30 +j2 

b.  0,  =  30”.  I  =  0.1  m,  o  =  0.2.  L  =  3T.  e  =  30  +  j2 

c.  01  =  30°,  X,  =  0.1  m.  o  =  0.5.  L  =  3T.  e  =  30  +  J2 


The  relative  roughness  of  a  surface  Is  a  luncUon  of  the  wavelength  of  the  Incident  wave  so 
that  the  frequency  of  the  Incident  wave  will  also  affect  the  statistical  distribution  of  the 
scattered  power.  Decreasing  the  wavelength  Increases  the  Rayleigh  roughness  parameter  Z,  so 
that  the  relative  heights  of  the  scatterers  Is  Increased.  Figures  17  through  20  show 
comparisons  for  cases  where  only  the  wavelength  Is  changed. 

In  Figure  17,  the  statistics  are  for  0i  =  30°.  a  =  0.2.  and  L  =  3T.  Figure  17a  Is  for  the  longer 
wavelength.  X  =  0.25  m  while  In  17b.  X  =  0.1  m.  Figure  18  Is  similar  to  Figure  17  but  for  a 
rougher  surface,  and  In  Figure  19  the  angle  of  Incidence  Is  Increased  to  75°.  Figure  20  Is  the 
same  as  Figure  19  but  lor  horizontal  polarization.  In  all  cases,  the  variance  of  the  scattered 
power  Is  seen  to  be  closer  to  Riiylelgh  when  the  Incident  field  has  a  wavelength  of  0.25  m.  For 
the  shorter  wavelength,  or  higher  frequency,  the  scattering  becomes  less  R<tyleigh.  Table  2 
shows  the  variance  of  the  scattered  field  for  several  wavelengths  normalized  to  the  Rayleigh 
value  for  Gj  =  30°,  o  =  0.2,  L  =  3T.  0g  =  0°.  and  Cr  =  30  +  J2-  Over  the  range  of  wavelengths  tested, 
the  scattering  becomes  less  Rayleigh  as  the  wavelength  decreases,  or  the  roughness  parameter 
of  the  surface  Increases. 


Table  2.  Variance  as  a  Function  of  Wavelength. 


X  (m) 

oo  ,  oo 

.25 

2.891 

3.772 

1.22 

.15 

2.939 

4.526 

1.44 

.10 

2.999 

5.038 

1.60 

.05 

2.987 

5.614 

1.83 
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Figure  1 7.  Mean  tind  Variance  of  Scattered  Power  for: 

a.  e,  =  30“,  \  =  0.25  m,  o  =  0.2.  L  =  3T,  E  =  30  +  J2 

b.  Oi  =  30“  A.  =  0.1  rn,  o  =  0.2,  L  =  3T.  e  =  30  +  j2 


>w(m)0.250 


e  18.  Mean  and  Variance  of  Scattered  Power  for: 


a(m)  0.500 
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Figure  19.  Mean  and  Variance  of  Scattered  Power  for: 

a.  0i  =  75°.  X  =  0.25  m.  o  =  0.5.  L  =  3T.  e  =  30  +  J2 

b.  e,  =  75°.  X  =  O.l  m,  o  =  0.5.  L  =  3T.  e  =  30  +  j2 


Figures  21  and  22  are  plots  of  the  normalized  variance  as  a  function  <jf  rms  surface  heights. 
These  values  are  calculated  by  dividing  the  standard  deviation  by  the  mean  value.  Because  the 
Rayleigh  variance  Is  defined  as  o°°  =  (<r°)^,  a  value  of  one  means  the  variance  Is  equal  to  the 
Rayleigh  variance,  while  a  value  greater  than  one  means  the  varlcince  is  larger  than  the 
Rayleigh  value.  These  plots  show  that  as  Gg  approaches  the  backscatter  region  the  variance  Is 
larger  for  the  smoother  surface.  This  Is  because  the  smoother  surface  will  have  so  few  facets 
with  slopes  large  enough  to  cause  backscatter  that  the  small  number  of  facets  becomes  the 
dominating  factor  in  determining  the  scattering  distribution  for  this  region.  In  comparing 
Figures  21  and  22  II  can  also  be  seen  that  the  variance  Is  larger  for  6i  =  7.5".  This  would  be 
expected  because  the  slopes  required  to  cause  backscatter  would  have  to  be  even  steeper  for  this 
case.  In  the  near  specular  and  forward  scatter  regions,  the  scattering  from  the  smoother 
surface  Is  closer  to  Rayleigh  than  the  rough  surface  scattering  because  there  are  many  more 
facets  oriented  favorably  to  cause  scattering  In  these  regions. 

Table  3  compares  the  calculated  mean  In  the  backscatter  direction  with  measured  values 
given  by  Ulaby^  and  Long.^  The  measured  backscatter  from  soil  given  by  Ulaby  is  compared  to 
the  calculated  value  for  a  dielectric  constant  of  Cr  =  30  +  J2.  The  measured  backscatter  from 
Arizona  desert  given  by  Long  Is  compared  to  the  backscatter  calculated  with  Cr  =  2  +  J1.6. 
Seawater  measurements  from  Long  are  compared  with  calculations  for  £r  =  80  +  J9. 

For  the  higher  sea  state  calculation,  corresponding  to  16-ft  waves,  tlie  shadowing  function 
was  Included.  For  pure  backscatter,  0i  =  0g,  the  shadowing  function  Is  given  as:® 


S  =  (l  +C2)-' 


where 


C 


2 


,  ,  ,1/2  (-T^  cot^  0 

( ^  (tan  0J  exp  I - ^ 


-  erfc 


T  cot  0^  'I 
2o  /■ 


(10) 


For  the  parameters  used  In  this  calculation,  one  obtains  S  =  0.928.  This  means  that  even  for 
high  sea  states  at  a  70"  Incident  angle,  shadowing  Is  not  a  significant  factor. 


^  Ulaby,  F.T..  Dol)Son,  M.C.  (1989)  Radar  Scattering  Statistics  for  Terrain.  Arlech  House. 
Dedham,  Massachusetts. 

®  Sancer,  M.I.  (1969)  Shadow  Corrected  Electromagnetic  Scattering  From  a  Randomly  Rough 
Surface.  IEEE  Trans.  Antennas  Propagat..  AP-17:577-585. 
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Figure  2 1 .  Comparison  of  Normalized  Variance  for  Different  Surface 
Roughness  for  Gj  =  75°,  A.  =  0. 1  m.  L  =  3T,  e  =  30  +  j2 


Figure  22.  Comparison  of  Normalized  Variance  for  Different  Surface 
Roughness  for  0;  =  30“.  X  =  0.25  m,  L  =  3T.  e  =  30  +  j2 


Table  3.  Comparison  of  Measured  o°  With  Calculated  a°. 


Er 

einc. 

o  (m) 

T  (m) 

L  (m) 

30° 

.1 

1 

3 

-15.96  dB 

-13.9  dB 

liBSfEEl 

30° 

.1 

1 

3 

-23.61  dB 

-25.75  dB 

BBWWI 

30° 

.61 

6 

12 

-19.94  dB 

-20.0  dB 

80  +  19 

70° 

1.225 

6 

18 

-21.82  dB 

-31.0  dB 

For  the  given  terrain  surfaces,  as  well  as  the  lower  sea  state,  the  agreement  between 
{ alculatlons  and  measurements  Is  extremely  good.  The  calculation  for  the  higher  sea  state, 
however,  does  not  agree  with  the  measured  value,  but  this  Is  not  unexpected.  High  sea  states  do 
not  have  a  Gaussian  height  distribution  or  a  Gaussian  correlation  function,  so  the  model  used 
In  this  study  cannot  accurately  describe  scattering  from  very  rough  seas. 


4.  ARBITRARY  LINEAR  POLARIZATION 


When  the  Incident  field  Is  linearly  polarized  at  some  arbitrary  angle,  \|i,  from  the 
horizontal,  the  calculation  of  the  reflection  coefficient  becomes  somewhat  more  complicated. 
The  Incident  linearly  polarized  wave  can  be  separated  Into  a  horizontally  polarized  component 
and  a  vertically  polarized  component.  The  horizontal  and  vertical  reflection  coefilclents  and 
the  magnitude  of  the  overall  reflection  coefficient.  |R^l|  can  then  be  calculated. 


r)  sln%  • 


(11) 


The  scattered  field  will  no  longer  be  linearly  pwlaiized  but  will  have  some  elllptlclty  due  to 
the  unequal  phase  change  of  the  horizontal  and  vertical  components.  The  polarization  of 
the  scattered  field  can  be  found  by  examining  the  complex  polarization  factor,  p.  defined  as: 


E  I  E 

P  =  --  =  cosq  +j 

I  I 


H 


Sin  (|) 


(12) 


where:  6  phase  difference  between  vertical  and  horizontal  components  of  the  scattered 

field 
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Re(p) 


cos  (> 


Im  (p) 


sin  <1). 


The  magnitude  and  phase  of  the  components  of  the  scattered  field  are  determined  directly  from 
Rv  and  Rh-  Table  4  shows  the  real  and  Imaginary  components  of  p  for  0i  =  30”,  A.  =  0. 1  m. 
a  =  0.2  m,  L  =  3T.  y  =  60°,  and  Cr  =  30  +  J2.  The  very  small  magnitude  of  Im  (p)  shows  that  the 
polarization  of  the  scattered  field  Is  very  close  to  linear.  Case  2  Is  for  =  ‘^  *  j '  6  with  y  =  50°. 
The  scattered  field  In  this  case  Is  slightly  more  elliptical  due  to  the  higher  relative  magnitude 
of  the  dielectric  loss  factor,  e",  but  It  Is  still  highly  linearly  polarized.  These  two  cases 
represent  the  highest  degree  of  elllptlclty  for  the  values  of  p  that  were  examined. 


Table  4.  Complex  Polarization  Factor 


0s 

Case  1 

Case  2 

Re  (p) 

Im  (p) 

Re  (p) 

Im  (p) 

-45 

-1.7212 

-0.0003 

-1.1678 

-0.0084 

-30 

-1.7331 

0 

-1.1918 

0 

-15 

-1.7212 

-0.0003 

-1.1678 

-0.0084 

0 

-1.6888 

-0.0015 

-1.0970 

-0.0326 

10 

-1.6550 

-0.0026 

-1.0255 

-0.0564 

20 

-1.6116 

-0.0039 

-0.9360 

-0.0849 

30 

-1.5582 

-0.0054 

-0.8304 

-0.1168 

40 

-1.6116 

-0.0039 

-0.9360 

-0.0849 

50 

-1.6550 

-0.0026 

-1.025!) 

-0.0564 

60 

-1.6888 

-0.0015 

-1.0970 

-0.0326 

75 

-1.7212 

-0.0003 

-1.1678 

-0.0084 

90 

-1.7321 

0 

-1.1918 

0 

Im  (p)  =  0 

Linear  polarization 

p  s:  oo 

Vertical  polarization 

Im  (p)  >  0 

Right  elliptical  polarization 

p  =  i 

Right  circular  poleU'izatlon 

Im  Ip)  <  0 

Left  elliptical  polarization 

p  =  -1 

Left  circular  polarization 

p  =  0 

Horizontal  polarization 

Case  1;  01  =  30°.  X  =  0.1  m,  o  =  0.2  m,  L  =  3T.  y  =  60°.  =  30  +  J2 
Case  2:  0|  =  30°.  X  =  0. 1  m,  o  =  0.2  m,  L  =  3T.  y  =  50°,  Er  =  2  +  J 1 .6 


34 


Figures  23  through  25  show  the  mean  and  variance  of  the  scattered  power  as  the  tilt  angle  \|/ 
of  the  linear  polarization  increases.  While  the  mean  does  decrease  as  the  polarization  goes 
from  horizontal,  y  =  0”,  through  several  values  of  y  and  finally  to  vertical  polarization,  = 
90°,  the  tilt  angle  of  the  polarization  does  not  affect  the  Rayleighness  of  the  statistics. 


5.  CONCLUSIONS 


The  main  ob|ertlve  of  this  study  was  to  determine  whether  cell  size  has  any  efl'ecl  on  the 
mean  normalized  cross  section  and  variance  of  the  scattered  power  from  a  rougli  surface. 
Several  other  variables  were  also  studied  Including  wavelength,  dielectric  constant,  angle  of 
Incidence,  surface  roughness,  and  polarization. 

Throughout  the  figures  shown,  several  general  trends  can  be  seen.  It  Is  quite  clear  that  the 
cell  size  does  affect  the  scattering  statistics.  For  the  largest  cell  size,  L  =  12T.  the  calculated 
value  of  o°°  is  almost  an  exact  match  to  the  Rayleigh  value  of  o°°  over  the  entire  range  of 
scattering  angles.  As  the  size  of  the  cell  decreases,  the  variance  of  the  scattered  power  becomes 
larger  than  the  Rayleigh  variance,  especially  In  the  backscatter  region,  due  to  fewer  facets 
scattering  In  this  direction.  This  means  that  the  problem  of  non-Rayleigh  scattering  for  finite 
clutter  cells  is  a  more  significant  problem  for  monostatic  radar  systems  than  for  large  angle 
blstatlc  radars. 

The  wavelength  also  affects  the  statistics  of  the  scattered  power  for  a  surface  containing 
many  scatterers.  As  the  wavelength  Increases,  the  Rayleigh  parameter,  Z.  decreases  and  the 
surface  becomes  less  rough,  so  that  the  scattering  is  closer  to  Rayleigh. 

Other  parameters  such  as  dielectric  constant  and  Ult  angle  of  the  Incident  polarization 
affect  the  mean  value  of  the  scattered  field,  but  do  not  affect  the  Rayleighness  of  the  scattered 
field. 
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3.  Mean  and  Variance  of  Scattered  Power  for: 


'Igure  24.  Mean  and  Variance  of  Scattered  Power  for: 
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